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Azimuthal anisotropy (1)2) and two-particle angular correlations of high pr charged hadrons have 
been measured in Au+Au collisions at ^/sjvjv = 130 GeV for transverse momenta up to 6 GeV/c, 
where hard processes are expected to contribute significantly. The two-particle angular correlations 
exhibit elliptic flow and a structure suggestive of fragmentation of high pr partons. The monotonic 
rise of V2{pr) for pr < 2 GeV/c is consistent with collective hydrodynamical flow calculations. At 
Pt > 3 GeV/c a saturation of V2 is observed which persists up to pr = 6 GeV/c. 

PACS numbers: 25.75.Ld 



Collisions of heavy nuclei at ultra-relativistic energies 
exhibit strong collective flow effects indicative of a vol- 
ume of hot matter so dense that descriptions involving 
hydrodynamic behavior in a locally thermalized system 
may apply Q. The azimuthal anisotropy of final state 
hadrons in non-central collisions t| is sensitive to the 
system evolution at early times 3J. At high px 1 a hy- 
drodynamic description of the system may break down 
as processes involving hard scattering of the initial-state 
partons are expected to play the dominant role. 

Calculations based on perturbative QCD predict that 
high energy partons traversing nuclear matter lose en- 
ergy through induced gluon radiation ,41] , where the mag- 
nitude of the energy loss is dependent upon the den- 
sity of the medium |5(. Recent measurements of inclu- 
sive charged hadron distributions in Au+Au collisions at 
v / s]v7V=130 GeV show a suppression of hadron yields at 
high pt in central collisions relative to peripheral col- 
lisions and scaled nucleon-nucleon interactions, consis- 
tent with the picture of partonic energy loss in a dense 
system 0, Q. The fragmentation products of partons 
that have propagated through the azimuthally asymmet- 
ric system generated by non-central collisions may ex- 
hibit azimuthal anisotropy due to energy loss and the 
azimuthal dependence of the path length, providing im- 
portant information about the initial conditions and dy- 
namics in a heavy ion collision [8|, |9j . 



The azimuthal anisotropy of an event in momentum 
space is quantified by the coefficients of the Fourier de- 
composition of the azimuthal particle distribution with 
respect to the reaction plane, with the second harmonic 
coefficient V2 referred to as elliptic flow. These coeffi- 
cients can be inferred from the particle distribution with 
respect to the estimated reaction plane orientation, cor- 
rected for the reaction plane resolution, or from two- 
particle correlation analysis The methods are iden- 
tical if the azimuthal correlation between particles re- 
sults solely from their correlation with the reaction plane. 
Correlations that are localized in both rapidity and az- 
imuthal angle are characteristic of high energy partons 
fragmenting into jets of hadrons. Such short-range cor- 
relations may be isolated from elliptic flow using two- 
particle correlation analyses performed in different re- 
gions of relative pseudorapidity. 

The transverse momentum dependence of i>2 has been 
previously measured in Au+Au collisions at y/SNN = 
f30 GeV for charged and identified particles in 
the region of pr < 2 GeV/c. Elliptic flow at RHfC can 
be described by a hydrodynamical model for px up to 
2 GeV/c. In this Letter we report the first results on 
w 2 (pt) of charged particles measured in this reaction up 
to pr = 6 GeV/c, together with the analysis of two- 
particle azimuthal correlations among high pt charged 
particles. 
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The Solenoidal Tracker at RHIC (STAR) consists of 
several detector subsystems in a large solenoidal mag- 
net. The main tracking detector is the Time Projection 
Chamber (TPC), which has wide acceptance in pseu- 
dorapidity and complete azimuthal coverage 0. For 
this analysis, the full data set from the first year data 
taking of the STAR experiment was used, consisting of 
300K minimum-bias and 400K centrally triggered events. 
The minimum-bias data contain hadronic Au+Au inter- 
actions at i/SATAr=130 GeV corresponding to ~90% of 
the geometric cross section a geo , while the centrally trig- 
gered data provide an unbiased event sample for the most 
central 10% of the minimum-bias data set. The trigger 
conditions and event and track selection cuts for v 2 anal- 
ysis are identical to those used previously |Tl| . 

The reaction plane analysis method involves the cal- 
culation of the orientation of the event plane, which is 
an experimental estimator of the true reaction plane an- 
gle. For this analysis, the second harmonic event plane 
angles, fy 2 , were calculated for the full event and two 
subevents, consisting of randomly selected particles from 
the same event. The results are insensitive to the selec- 
tion method (random, pseudorapidity, or charge sign) for 
assigning particles to subevents Q. 

Jets may bias the reconstruction of the reaction plane 
if the intra-jet correlations produce asymmetries which 
are of similar magnitude to that due to collective flow. 
Systematic studies were undertaken to assess the bias. 
In general, the products of jet fragmentation have higher 
Pt than other particles produced in a collision |13| . Only 
low pt particles (pr < 2 GeV/c) were selected to cal- 
culate the event plane orientation. Using p-p-cutoffs of 
1.5, 1.0 and 0.5 GeV/c resulted in statistically consistent 
v 2 values. In addition, two sets of events were analyzed, 
the first containing only events with a particle of pt > 3 
GeV/c and the second with all events, separately. The 
values of V2 as a function of pr were the same in both data 
sets. To ensure that particles produced within a jet do 
not affect the reaction plane reconstruction, all subevent 
particles in a pseudorapidity region of |A?7| < 0.5 around 
the highest pr particle in the event were excluded and 
the results were insensitive to this procedure. Varying 
the track selection criteria (distance of closest approach 
to the primary vertex, number of measured space points, 
etc.) also made no significant difference. 

Finite momentum resolution at high transverse mo- 
menta combined with a rapid decrease of hadron yield 
with increasing px may cause flattening of the pt- 
dependence of v 2 . The momentum resolution has been 
determined by embedding simulated single tracks into 
real raw data events. For the cuts used in this analysis, 
the momentum resolution for px > 1.5 GeV/c is parame- 
terized as Sp T /p T = 0.013 + 0.014p T /(GeV/c) for central 
events at magnetic field _B=0.25 T. We have studied the 
possible effects of the momentum resolution on v 2 at high 
transverse momenta using a Monte-Carlo simulation by 



generating particles with a power-law py-distribution Q 
and with various V2{pr) dependencies. The estimated 
relative systematic error on v 2 due to momentum resolu- 
tion is 5% at pt—5 GeV/c. 

The reaction plane analysis integrates all possible 
sources of azimuthal correlations, including those un- 
related to the orientation of the reaction plane. The 
non-flow correlations may be due to resonance decays, 
(mini)jets, final state interactions (particularly Coulomb 
effects), momentum conservation, etc. The strength 
of non-flow contributions relative to the measured az- 
imuthal asymmetry was discussed in Q] , which estimated 
that 15-20% of the v 2 signal obtained with a reaction 
plane analysis method is due to non-flow correlations. A 
four-particle correlation method |l4| for flow measure- 
ments reduces non-flow sources to a negligible level. The 
centrality averaged values of v 2 from four-particle cor- 
relations are 15% lower than those obtained from the 
reaction plane analysis [l5j |. The ratio of v 2 from the 
two methods is approximately independent of pseudora- 
pidity and transverse momentum within 0.1 < px < 4.0 
GeV/c, the range accessible to the four-particle correla- 
tion method with the current statistics [l^. Based on 
these studies, we assign a ^~ 20 % systematic uncertainty 
to the v 2 values presented. This uncertainty is highly 
correlated point-to-point at all pr-values. 

Figure ^ shows the azimuthal distributions with re- 
spect to the reaction plane of charged particles within 
2 < pt < 6 GeV/c, for three collision centralities. The 
distributions are corrected for the reaction plane reso- 
lution. All distributions exhibit the second harmonic 
behavior characteristic of elliptic flow. A fit of dis- 
tributions by 1 + 2v 2 cos2((j) lab - * p ; ane ) yields v 2 = 
0.218 ± 0.003, 0.162 ± 0.002, 0.090 ± 0.001 for the three 
centrality bins. The errors are statistical only. There is 
large azimuthal anisotropy at higher for all centralities. 

Figure [21 shows the differential elliptic flow v 2 as a 
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FIG. 1: Azimuthal distributions with respect to the reaction 
plane of charged particles within 2 < pr < 6 GeV/c, for three 
collision centralities. The percentages are given with respect 
to the geometrical cross section o geo . Solid lines show fits by 

1 + 2v 2 COs2(0 i(lb - 'S! plane)- 
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FIG. 2: V2(pr) for different collision centralities. The errors 
are statistical only. The systematic uncertainties, which are 
highly correlated point-to-point, are ^^Vo. 

function of px for three collision centralities. The val- 
ues of V2(pt) were obtained from the second moment of 
the distribution of particles with respect to the reaction 
plane, i.e. from the average values (cos2(0; a b — ^piane)}, 
corrected for the reaction plane resolution. At a given 
Pt, the more peripheral collisions have larger v^. For 
all centralities, v-i rises linearly up to px = 1 GeV/c, 
then deviates from a linear rise and saturates for px > 3 
GeV/c. The saturation persists up to 6 GeV/c and is 
in contrast to non-dissipative hydrodynamical calcula- 
tions, which predict a continuous rise of t>2 with increas- 
ing transverse momentum |l6| . 

In Ref. [9j , the particle production is decomposed into 
phenomenological "soft" and perturbative QCD calcu- 
lable hard components. The soft nonperturbative com- 
ponent incorporates hydrodynamic elliptic flow, whereas 
the pQCD calculable part includes energy loss (jet 
quenching). In this model the magnitude of t>2 at high 
Pt is sensitive to nuclear geometry and the initial gluon 
density achieved in a collision. Figure [3] compares the 
minimum-bias differential elliptic flow V2(pt) with calcu- 
lations from Ref. Q. These calculations also predict a 
decrease of V2 with increasing pt at high transverse mo- 
menta. This decrease should even be stronger if trans- 
verse expansion of the system is taken into account |l7| . 
A rapid expansion dilutes the initial coordinate space az- 
imuthal asymmetry resulting in a reduction of the mea- 
sured azimuthal anisotropies due to energy loss. A flavor 
dependence of both the asymmetry and pr differential 
particle multiplicities has been suggested as one of the 
possible scenarios for the V2{pt) behavior 01- The ob- 
served saturation of V2 at p? ~2-3 GeV/c can be quanti- 
tatively reproduced in a parton cascade model with only 
elastic rescatterings, but extreme initial gluon densities, 
dN 9 jdr\ ~ 15,000, or extreme elastic parton cross sec- 
tions, ~ 45 mb, are required fljj ]. 

In order to verify the existence of a hard scattering 
and fragmentation component at high pt, two-particle 
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FIG. 3: V2(pr) for minimum-bias events (circles). The er- 
ror bars represent the statistical errors and the caps show 
the systematic uncertainty. The data are compared with hy- 
dro+pQCD calculations |3] assuming the initial gluon den- 
sity dN 9 /dy = 1000 (dashed line), 500 (dotted line), and 200 
(dashed-dotted line). Also shown are pure hydrodynamical 
calculations (solid line). 

angular correlation measurements are used [liij . In cen- 
tral heavy-ion collisions, it is impossible to reconstruct 
jets fully due to the large overall particle density. An- 
gular correlations of high pt particles, however, allow 
for an identification of a hard scattering component on a 
statistical basis. The fragmentation of high pt partons 
into several particles results in correlations of hadrons 
at small At), A<j>. In order to isolate this short-range 
component of the two-particle correlation function, the 
azimuthal correlations of high px particles are measured 
in two regions of relative pseudorapidity. At large A77, 
we assume that the azimuthal correlations are free from 
the fragmentation component. 

For the two-particle azimuthal correlation analysis, 
events containing a trigger particle having 4 < p^(trig)< 
6 GeV/c and \r)\ < 0.7 are used. For these events, 
we measure the relative azimuthal distribution of other 
charged tracks with 2 GeV/c < pt < Pt (trig) and 
1 77 1 < 0.7, and the distribution is normalized to the num- 
ber of high pt trigger particles, 

N 1 = tt~~ ~ / dAriN(A4>, A 77). (1) 

^trigger d[A<p) iVtrigger C J 

■^trigger is the observed number of tracks satisfying the 
trigger requirement, and N(A<f>, Arj) is the number of 
observed pairs as a function of relative azimuth (A</>) 
and pseudorapidity (A77), and e is the single track effi- 
ciency. Due to the nearly uniform azimuthal acceptance 
of STAR, no mixed event reference is required. With 
this definition of the correlation function, the efficiency 
for finding the trigger particle cancels, and we need only 
correct the data for the efficiency of finding the lower pt 
particle. This efficiency is determined using embedding 
and is 66% for the tracks used in this analysis. 
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FIG. 4: High pr azimuthal correlation functions for central 
events. Upper panel: Correlation function for |Ary| < 0.5 
(solid circles) and scaled correlation function for 0.5 < | Arj\ < 
1.4 (open squares). Lower panel: Difference of the two corre- 
lation functions. Also shown are the fits to the data (described 
in the text). 

Figure 2| shows the two-particle azimuthal correlation 
functions for |A?7| < 0.5 and |A?7| > 0.5 for central 
events. Here the small relative pseudorapidity correlation 
function is absolutely normalized, while the large rela- 
tive pseudorapidity correlation function has been scaled 
to match the small \Arj\ correlation function in the re- 
gion 0.75 < |A0| < 2.25. There is an enhancement near 
A(f> = 0. This short range correlation may be evidence 
of hard scattering and fragmentation. The large relative 
pseudorapidity correlation function is fit by a functional 
form dN/d(A<fi) oc 1 + 2v\ cos(2A</>), expected if the cor- 
relations are due entirely to elliptic flow. This fit gives 
v 2 = 0.11 ±0.02. 

In order to compare this pair-wise method quantita- 
tively to the reaction plane-based v 2 analysis shown in 
Figure ^ we construct azimuthal correlation functions 
with a trigger particle having 2 < pr < 6 GeV/c. 
For |A?7| > 0.5, the azimuthal correlations give v 2 re- 
sults consistent with the reaction plane analysis: v 2 = 
0.203 ± 0.012,0.160 ± 0.007, and 0.091 ± 0.003 for the 
three centrality bins. If no pseudorapidity gap is required 
between two particles, the v 2 values are 10-15% higher. 

Also shown in Figure 0] is the difference between the 
small and large relative pseudorapidity correlation func- 
tions. There is an enhancement near zero, and a flat 
correlation function at larger A<j>. A Gaussian fit gives 
a = 0.27 ± 0.09(stat.)±0.04(sys.) radians, where the 
systematic error was estimated by varying the binning 
and range used to scale the correlation function for 
0.5 < \Ar)\ < 1.4. This width is consistent with ear- 
lier observations of jet characteristics in pp collisions at 
slightly lower 0] and higher 0] energies. The HIJING 
event generator [2(|, where hard scattering and frag- 
mentation dominate particle production at these trans- 
verse momenta, predicts a = 0.20 ± 0.01 for the same 



kinematic cuts.. Integrating the signal we observe that 
4.9 ± 1.7(stat.)±0.4(sys.) % of charged particles with pr 
of 4-6 GeV/c have an associated charged particle with 
Pt > 2 GeV/c. The systematic error is dominated by the 
uncertainty in the absolute efficiency determined via the 
embedding procedure. The contributions of resonance 
decays and photon conversions were studied and found 
to be insignificant. 

In summary, the measurements of azimuthal anisotro- 
py v 2 of charged particles with pt of 3-6 GeV/c reveal 
a saturation pattern of v 2 with values that decrease sys- 
tematically with increasing centrality. This contradicts 
non-dissipative hydrodynamics which predicts a mono- 
tonically increasing v 2 with increasing px, but the data 
may be consistent with dissipative dynamics with finite 
parton energy loss. In addition, a comparison of the two- 
particle azimuthal correlation functions for particles with 
|A?7| < 0.5 and |A?7| > 0.5 suggests the existence of a 
short-range correlated component at high px in addition 
to underlying global elliptic flow. This may be the first 
direct evidence at RHIC for hard scattering and parton 
fragmentation. The data provide important constraints 
on the theoretical interpretations of the mechanism of 
high pt particle production in ultra-relativistic heavy- 
ion collisions. 
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